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Abstract 



We propose a method which allows the generalization of the Landau lattice 
gauge-fixing procedure to generic covariant gauges. We report preliminary nu- 
merical results showing how the procedure works for SU{2) and SU{2>). We also 
report numerical results showing that the contribution of finite lattice-spacing ef- 
fects and/or spurious copies are relevant in the lattice gauge- fixing procedure. 



1 Introduction 



Recently, lattice QCD Monte Carlo simulations requiring a gauge-fixing have be- 
come relevant; the gauge-fixing is essential to study gauge dependent quantities 
like, for example, the propagators of the fundamental fields entering the contin- 
uum QCD lagrangian and to use smearing techniques Q. 

The study of gluon and quark propagators, among other things, allows a bet- 
ter understanding of the infrared behaviour of the theory and the confinement 
mechanism. Moreover, quark and gluon matrix elements can be used to obtain 
renormalization conditions as proposed in |§, ^ . 

In the last few years, numerical studies of lattice propagators have been performed 
by several groups: the aim of the authors was to study the mechanism 

through which the gluon may become massive at long distances whereas more re- 
cent attempts studied its behaviour as a function of momentum 0]-^]. Analogous 
studies of the quark propagator exist 

The existence of lattice Landau and Coulomb gauge-fixing ambiguities has been 
verified |10]-[^]. Studying the characteristics of these ambiguities and their influ- 
ences on gauge fixed quantities is interesting at least for two reasons: the existence 
of these ambiguities could be the analogous of the Gribov problem in the contin- 
uum formulation of non abelian gauge theories |14] and the gauge-fixing is essential 
for the analytical study of the continuum limit of lattice gauge theories. 
In practice, there are some cases in which it is convenient to use a gauge de- 
pendent procedure to compute gauge invariant quantities. For example, smeared 
fermionic interpolating operators are being used in lattice QCD spectroscopy and 
phenomenology. The smearing operators are gauge dependent and therefore the 
gauge must be fixed before they are calculated. 

Up to now, the only covariant gauge for which it is known the algorithm to fix 
it on the lattice is the Landau Gauge |15|-[^]. This algorithm uses the original 
idea of Gribov |14| in the continuum, restricting the domain of integration of the 
partition function in the region where the functional F{Q) = reaches an 

extreme. Also algorithms for non covariant gauges are discussed in literature (see 
for example p^). 

In this paper we propose a procedure which allows to generalize the Landau gauge- 
fixing on the lattice for a generic covariant gauge. This gauge-fixing procedure 
could verify the gauge independence of some results obtained in literature |Q]-||9| 
and could allow to discriminate between gauge artifacts and true physical proper- 
ties of the fundamental fields entering the QCD lagrangian. 
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2 Covariant gauges quantization 



In this section we will briefly review the general formalism to quantize a non abelian 
gauge theory using covariant gauge conditions ]l9[ . 

Neglecting the Gribov problem, let us assume that we can find a gauge section 
in the space of gauge fields which intersects once and only once all gauge orbits. 
The Landau gauge is readily extended to an auxiliary gauge-fixing condition of the 
form 

d,A^^{x) = A{x) (1) 

where A(x) belongs to the Lie algebra of the group. Since gauge-invariant quanti- 
ties should not be sensitive to changes of gauge condition, it is possible to average 
over A(x) with a gaussian weight 

ZiJo)= JsAe-^I''"^^'^''"^ J M^<57?5f/e-^(^)-Vo.*('7,^.A.)+/^oCP5(9^^^_A) 

(2) 

obtaining the standard formula 




where O is a gauge-invariant operator. 

The choice a = 1 is referred to as Feynman gauge instead for q = the Landau 
gauge is recovered. 

In the next sections we will show how to implement such a formulation for a non 
perturbative numerical simulation to compute the mean value of a gauge-dependent 
observable on the lattice using the expression (|2|) for Z{Ja). 

3 The functional for covariant gauges 

The functional proposed by Gribov, directly in the continuum, in order to fix the 
Landau Gauge is 

F{n) ^\\A^\\' = Jtt {A^{x)Ai{x)) d'x (4) 

where 

A^{x) = VL{x)A^{x)Vt^{x) - Ul{x)d^,^\x) . (5) 

and 

^{x) = e^"'^^) (6) 
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is a group matrix, wi^x) and A^{x) belonging to the Lie algebra of the group. 
In appendix A we show that 

^ = --^id,A^,r^^\v.) (7) 



with 



Equation shows that F{Q) is stationary when d^A^ = 0. 

In order to fix the gauge discussed in the section ^ we should be able to find a 
functional H{0,) stationary when 

d.A^ix) = Aix) , (9) 

with A(x) having a gaussian distribution. The most naive way to define H{0,) 
would be to find a functional h{0,) such that 

= ^ {a-{xW\w{x))) (10) 



so that 

= -| {d.A^^ix) - A{x)r'^'^\w{x) 

However we will now show that this is not possible; in fact for a non abelian gauge 
theory does not exist a functional satisfying (p!o|). A necessary condition for the 
existence of such a functional would be 



5w'^(x)5w^{y) 5w^{y)5w^{x) 
which implies the integrability condition 



^ -{A''{y)^''\w)) = — J--(A<^(x)$°"(u;)) . (13) 



5w'^{x) 6w^{y) 

Expanding ^"'^(w{x)) in power of 'w{x), equation ( p^ ) should be satisfied order by 
order in w{x). From equation (P) we have 

<^''\w) ~ 5"'' + I"'' = 5"^ + Y , (14) 
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Equation (|T^) is then in contrast with the antisymmetry of f°'^'^. 

The functional we propose in order to fix the gauge discussed in section |2| is: 



In fact 



5H{n) 



{d.A'l - A){d,A^ - A) 



5w°- 



which shows that H{Q) is stationary when 



aX(^)-A(x) = o . 



(15) 
(16) 

(17) 



The problem is to show that no spurious stationary points of ( p]5| ) exist. In fact, 
for example, apparently 



df,A^/^{x) - A(x) = cost. 



also satisfies 



5H{VL) 

5w°- 







(18) 



(19) 



and other non trivial solutions of the equation ([T9|) could exist as well. In our 
case we can exclude the possibility of constant zero modes. In fact integrating the 
equation (^) we obtain 

J df,A^{x)d'^x - J k{x)(fx = cost J (fx . 

Through the use of periodic boundary conditions we have 

A{x)d'^x = cost J dt^x . (20) 

Equation (pO|) assures that cost = when A(a;) has a gaussian distribution. This 
result assures that the functional H{^) does not have spurious stationary points 
satisfying equation ([l8|). 

It is interesting to note that for A{x) = the functional H{Q) has absolute minima 
for any Q satisfying df^A^{x) = 0. This is not the case for the Gribov functional. 



4 Covariant gauges on the Lattice 

The gauge variables of a compact lattice gauge theory are the links U^{x) and they 
are elements of the gauge group. 

In this section we outline a procedure to compute numerically the mean value 
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of a gauge dependent operator O on the lattice using a generic covariant gauge 
quantization. The mean value of a gauge dependent operator on the lattice is 
defined as : 

{O) = y"(5Ae-^/"''^^''(^') J dUAf{U)6{A - A)e-^'^'^^^0{U) . (21) 

where SwiU) is the Wilson lattice gauge invariant action, Aj(C/) is the Faddeev- 
Popov determinant and 

I 4 

= 17- - Ul{x)]traceless " [U^i{x - /i) - C/^(x - lj)]traceless ■ (22) 



The computation of the integral (|2l| ) can be schematized as follows: 

• A gauge configuration {U} with periodic boundary conditions according to 
the gauge invariant weight e~^^^^^ is generated; 

• For each {[/} configuration random matrices A(x) belonging to the group 
algebra are extracted according to the weight e~2^'^^^ 

• Given A(x), a numerical algorithm extremizes a discretization of the func- 
tional H{^D,). This defines the lattice gauge-fixing condition 

A(x) - A(x) = Vx ; (23) 

• the mean value of the lattice gauge dependent operator is then defined as 

(O)^'^" = ^ E . (24) 

{conf} 

Since A{x) has a gaussian distribution for all x, if the number of lattice sites is 
sufficiently large 

i^A(x)=A = 0. 

X 

5 Lattice gauge-fixing for covariant gauges 

A possible discretization of H(Q) is 

i/i(J7) = iTV^[A^(x)-A(x)]' (25) 

X 

where 

1 4 

A^(X) = -—Y.[Uf}{x)-U^Hx)]traceless-[U^{x-fi)-U^Hx-fl)]traceless , (26) 
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A(x) are matrices which belong to the group algebra, V is the lattice volume and 
links transform under a gauge transformation as 

u^ix) = n{x)u^{x)nHx + fi) . 

H^iP) is the simplest discretization of H{0,). For a generic transformation close 
to we define 



6Q 



and we have 



K\x) + Uf\x - + Uf{x - ^)) + 



tQ'- 

V ' ""/i 



.(A^\x - ^) - A(x - M))(f/f - ^) + Uf{x - ^,)) + (27) 



(A^'^x + /i) - A(x + /.))(< t(^) + ufix)) 
where T"" are the group generators. 

We have to find a gauge transformation Q such that {Uj}} satisfies 

6Hl 



(5e'^(x) 



V X . (28) 



In order to solve numerically equation (^) we use an iterative algorithm for the 
minimization of Hl^Q). If such an algorithm converges (this is not guaranteed and 
must be checked in practice), its fixed points will be configurations satisfying the 
condition (|2|). 

In order to study the convergence of the algorithm, two lattice quantities can be 
monitored 

Hl = :^TY^[A(x)-A(x)f 

The function Hi{Vt) is defined on a compact set. If the numerical algorithm of 
minimization converges, 'dn must go to a value of the order of the precision required 
for the minimization whatever starting configuration we use. The procedure works 
if also Hl goes to a value of the order of the precision required for the minimization 
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for those configurations where the algorithm converges. In the next section we will 
show numerical results suggesting that Hl{Q) goes to zero when it is minimized 
for the configurations considered. 

We remark that the same discretized definition of gauge field {l/2iag)[Ui^i{x) — 
Uji{x)]traceiess is uscd in the definition of Hl{Q) and of A^(x). 
On the contrary in the Landau lattice gauge-fixing procedure, usually adopted 
in the literature, different definitions of the lattice gauge potential are used . In 
section ^ we report numerical results showing interesting phenomena arising if 
different lattice gauge potential definitions are used in the lattice Landau gauge- 
fixing procedure. 



6 Numerical Simulations 

In this section we report preliminary results of some numerical simulations. In this 
exploratory study we considered SU{2) and SU{3) non thermalized configurations 
with links defined as 

with 



where B^{x + ^) are smooth functions of x (typically sin(x) and cos(x)), T" are 
group generators and "6" is a parameter which determines the distance of U^{x) 
from identity. We generated, on a volume V = 4"^, nine configurations for each 
group, three for each value b = 0.1, b = 0.3, b = 0.7. Since the function Hl{Q) 
is defined on a compact set, the procedure works whatever configuration we start 
from, except for numerical convergence problems. A complete study of this nu- 
merical algorithm for thermalized Monte Carlo configurations will be presented in 
a future paper. 

The gauge-fixing algorithm implements an iterative minimization for Hl{0,) which 



updates link matrices via SU{2) subgroups, as proposed in ref.[20|. We monitored 
the quantities Hi and "dn after every lattice sweep of the gauge-fixing algorithm. 
We mention that after each gauge-fixing sweep we orthogonalized the link vari- 
ables, to make sure that they have not been driven off the gauge group. 

In fig. |l| we plot Hl and 'dn as a function of the number of gauge-fixing sweeps 
for a configuration of SU{3) with a = 0.3 and A{x) = 0. It is easy to see that 
both Hl and "dn go to values consistent with the precision required for gauge- 
fixing (< 10~^). For all configurations we generated the behaviour of and 
is consistent with that shown in fig.||. For A(x) = we fixed the gauge condition 
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Figure 1: Values of Hl and as functions of number of iterations Nuer of the covariant 
gauge-fixing algorithm obtained extremizing Hi for a configuration of SU{3) with a = 
0.3, y = 4^ and A{x) = 0. 

A(x) < 10~^ for each configuration also using the old Landau gauge- fixing algo- 
rithm. For the two configurations obtained with the two algorithms, we measured 
the discretization of -F(O) 

FLm = -^J2^eTrUl}{x) (29) 

x,ii 

as a measure of the distance between the two configurations. The values are the 
same with a precision of 10~^, consistent with the precision required to minimize 
Hl . 

In fig. ^ we plot Hl as a function of the number of gauge-fixing sweeps for a 
configuration of SU{2) with a = 0.1 and A(x) ^ 0. 
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Figure 2: Values of Hl and 'du as functions of number of iterations Nuer of the covariant 
gauge-fixing algorithm obtained extremizing Hl for a configuration of SU{2) with a = 
0.1, y = 4^ and A{x) ^ 0. 

In fig. § it is easy to see that the behaviour of the algorithm when A(x) 7^ is 
similar to that one with A(a;) = . 

Finally we stress that the numerical minimization of the functional ( p5|) is more 
complicated than that of Fl{^1)- In this exploratory study the time per iteration 
required to minimize Hl{^) is 5 — 10 times the time per iteration required to 
minimize Fl{^)- 
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7 Comments on the lattice Landau gauge-fixing 



In this section we report some numerical results which show some consequences 
of using different lattice definition of gauge potential in the Landau lattice gauge- 
fixing procedure. 

The standard way of fixing the Landau gauge on the lattice is based on the mini- 
mization of the function Fl{Q,) defined in the equation ( p9|) which is a particular 
discretization of F{Q) with 

IK]''^' = - 7^ (U, + Ul- 21). (30) 
g a 

Once Fl{^}) has been numerically minimized, the quantity 

^^ = ^ETr[A(x)At(x)] (31) 



X 



is supposed to vanish and this is the signal that the lattice gauge condition A(x) = 
has been achieved at each site. A(x) is a discretization of df^A^{x) with 

traceless 

^^^^ = 2i^g ■ ^"^^^ 

Some authors showed that, in the lattice Landau Gauge, lattice copies characterized 
by different values of exist. These solutions of A^(x) = are indeed different 
gauge-related configurations on the lattice not connected by a colour rotation. 
However it is important to note that the definition of the lattice gauge potential 
used in the discretization of F{Vl) does not correspond to the definition used in 
the discretization of the divergence A(x). This implies that two different lattice 
definitions of have been used to check the Landau gauge-fixing condition and 
to tag two different solutions. 

The authors in ref. ||2^, ^ have found two ensembles of " lattice copies" with V = 
16^ X 32 and /3 = 6.0. To analyse these configurations we used the usual Fl and 
dp variables defined in equations (^J3^) and another function defined as 

^F = :^E^[^'(^)^'^(^)] (33) 

X 

where 

1 4 

^'(^) = J—Y.\iU^.{x)-Ul{x)){U^{x)+Ul{x))]traceless 

- [{U^{x - H) - Ul{x - ll)){U^{x - H) + Ulix - ^l))]traceless 
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We stress that A and A' are two discretized definitions of d^Af^ wliich tend to tlie 
same expression as a ^ 0. For each "lattice copy" we measured the values of Fl, 
"dp and "d'p which we report on table |l[ These numerical results show that the 



Ensemble 


Copy 




dp 




A 


1 


2.583760909874 


2.006239333971^ 


- 10 


0.4114643051772 




2 


2.584445214939 


3.0326250142058£; 


- 10 


0.4090001249674 




3 


2.584426547617 


8.4843522906603£; 


- 10 


0.4091814689158 




4 


2.583991349769 


8.2208712589082£; 


- 10 


0.4106076010266 




5 


2.584036530091 


2.8552454724294^; 


- 10 


0.4097852320118 




6 


2.58354405348 


3.3687550405651S 


- 10 


0.4117689556597 


B 


1 


2.581884641098 


2.5513621387496^; 


- 10 


0.4177011836703 




2 


2.582231848321 


7.5376561548352^; 


- 10 


0.4159162059697 




3 


2.58231860688 


6.1085384335158S 


- 10 


0.4159582027108 



Table 1 : Final values of dp and d'p after gauge- fixing algorithm which extremizes numer- 
ically the functional for two gauge fixed ensemble of configuration with V = 16^ x 32 
and P = 6.0. 

finite lattice spacing effects and/or spurious copy contributions to "i^^ and i^'p are 
of the order of 10~^ while the difference between the values of Fi for two different 
copy are of order 10~^ when < 10"^''. This shows that with this method it 
is not possible to decide if two different solutions of Landau lattice gauge-fixing 
condition correspond to different Gribov copies in the continuum. This procedure 
is not apt to decide if the lattice multiple solution problem of A^(x) = has an 
analogy with the continuum Gribov problem. 

Moreover we observe that for a lattice copy there is a big difference between the 
values of the two discretizations i!}f and 'd'p respect to the precision required for 
the minimization. It would be interesting to understand if the difference between 
'dp and "dp is only due to the higher order lattice spacing effects or to spurious 



copies [25]. The problem of higher order contributions and/or spurious solutions to 
i^p afflicts all numerical computations of gauge dependent operators. In matching 
numerical results obtained on the lattice with the corresponding continuum for- 
mulas one must carefully evaluate the error assigned to the gauge-fixing condition 
even if < 10"^". Moreover the residual gauge freedom associated to lattice copies 
tagged with the functional ( p9|) should not induce effects higher than the sistematic 
uncertainty due to the higher order contributions and/or spurious solutions. 
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The numerical results shown in this section indicate also the importance to improve 
the gauge-fixing algorithm on the lattice also for the Landau gauge. 

Conclusions 

We have proposed a method which allows the generalization of the Landau lattice 
gauge-fixing procedure to generic covariant gauges. We have shown that a func- 
tional whose stationary points are dfj,A^{x) = A(x) cannot be obtained as a direct 
generalization of the Landau lattice gauge-fixing functional used in literature. In 
the continuum we proposed a functional reaching an extreme when d^A^{x) = A(x) 
and we used the simplest discretization of H{^}) to fix numerically a generic co- 
variant gauge on the lattice. We reported preliminary numerical results showing 
how this procedure works for SU{2) and SU{3). Numerical results also show that 
the contribution of finite lattice-spacing eff'ects and/or spurious copies are relevant 
in the lattice gauge-fixing procedure and must be carefully evaluate. 
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Appendix A 



In this appendix we calculate the variation of a SU{N) matrix 



(34) 



w 




a 



and of the functional 




(35) 



for an infinitesimal variation of w [^]. Following Feynman |22], ^}{w + dw) is 



0,(11) + 5w) 




12 



+ 1 / e 

10 

»1 

Q{w) + i I ds' 







1 s' 

i I . dsws r i I dsws 



Vl{w) + i ds' 
Jo 

Jt + iJ^ dse-'''"T^e'''"5w'' ) n{w) 



If we define 



and (t"*(s) is such a way that 

then it is obvious that 
and then 

If we define 7''^ = P^'^w'' then 
and then 



= -i[u;,S"(s)] 



We can conclude that 



el -1 



-I ab 



As is a unitary matrix 



(36) 
(37) 

(38) 

(39) 
(40) 

(41) 
(42) 

(43) 



(44) 



substituting the ( ^3| ) in the last equation and considering only linear terms in 5w 
we obtain 

[^^^\w)tHw'')^ = ^''\w)tHw'' (45) 
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and then 



(46) 



If we remember that 




n{x)Af,{x)n^x) - -n{x)di^n''{x) 



(47) 



n{x) 



(48) 



it is easy to verify that 



SF{n) = -- [ d^x(a^A^)«$"^(i(;)(5w;^ 



(49) 



and then 



dF{Q) 



(50) 
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